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Abstract 

We investigate the dynamics of two-photon correlations generated by the interaction of a three-level atom in 
the H, A or V configuration, with two classical external driving fields, under the rotating-wave approximation, 
in the presence of level decays. Using the example of a rubidium atom in each configuration, with field 
strengths validating the single-photon approximation, we compute measurement based correlations, such 
as measurement induced disturbance (MID), quantum discord (QD), and quantum work deficit (WD), 
and compare the results with that of quantum entanglement (concurrence). Certain correlation properties 
observed are generic, model independent and consistent with known results, e.g., MID is an upper bound 
on QD, QD and WD are monotonic, and the generic correlation behavior is strongly affected by the purity 
of the photon states. We observe that the qualitative hierarchy, monotonicity and steady-state behavior of 
the correlations can be controlled by the choice of parameters such as atomic decay constants and external 
driving field strengths. We point out how particular configurations are better suited at generating monotonic 
correlations in specific regimes and how the steady-state correlation behavior and hierarchy are affected by 
the population dynamics of the density matrix for different parameters. The possibility of using well studied 
quantum optical systems such as the three-level atom to generate, characterize and parametrically control 
mixed state quantum correlations establishes an important step in the direction of their implementation in 
quantum information tasks. 
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1. Introduction 

The interaction of atomic systems and external electromagnetic fields is a principal source for the gen- 
eration and classification of quantum correlations [1]. The quantum nature of these atom-photon systems 
and the ability to implement such systems in controlled experimental settings make them important tools in 
the study of nonclassical features [2]. From the perspective of quantum information theory (QIT), atomic 
systems are the quintessential computational hardware needed for the future implementation of quantum 
information protocols [3, 4], and photons are the basic building blocks of quantum communication [5] and 
cryptography [6], Hence, the generation and manipulation of nonclassical correlations in complex atomic 
systems interacting with radiation fields is one of the most challenging aspects of future applications of QIT. 

The various popular indicators of nonclassicality and measures of quantum correlations are dependent on 
the theoretical perspectives invoked for their quantifications and are often not consistent when the interacting 
states have sufficient amount of mixedness [7] . The main dichotomy in the definition of quantum correlations 
arises from the question of what constitutes quantumness. The extensively studied entanglement-separability 
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Figure 1: (Color online) A three-level atom in the (a) S, (b) A, and (c) V configuration. Ti and T2 are the decay constants 
(numbers shown in units of MHz) of the levels \2) and |3). ui, i>2, and f2i, O2 are the optical frequencies and the Rabi 
frequencies of the two near-resonant driving fields, oji and UJ2 are the two atomic transition frequencies. Ai and A2 are the 
field detunings, set to zero in this work. 



[8] criterion to define quantum correlations stems from the understanding that the main feature of quantum- 
ness arises from the superposition principle [9]. Another feature to have received widespread attention in 
recent times is the definition of nonclassicality or quantumncss on the basis of measurement based correla- 
tions. Such definitions take into account the fact that an important feature of quantumness in nature arises 
from noncommutativity of operators [10]. Any physical measurement on a quantum system disturbs the 
noncommutative nature of the system and thus effectively erases quantum correlations. Information theo- 
retic measures such as quantum discord (QD) [11, 12], quantum work deficit (WD) [13] and measurement 
induced disturbance (MID) [10] are based on the unique role of measurement in quantum physics. 

The nonclassical properties of three-level atomic systems have been well studied in quantum optics for 
understanding quantum-coherence phenomena such as electromagnetically-induced transparency (EIT) [14], 
lasing without inversion [15], and coherent trapping [16]. Three- level atoms interacting with low-strength 
driving fields, similar to EIT systems, have been used to generate entangled two-mode photon states which 
can be suitably manipulated to yield desired correlations [17]. The knowledge of nonclassical correlations 
carried by emitted photons in atomic systems may prove immensely useful in designing future QIT systems 
for communications and computation. The preparation and characterization of mixed quantum correlated 
states is a very challenging prospect with regards to generation of quantum states. Well known processes 
like parametric down-conversion [18] generate pure entangled states with poor conversion efficiency and few 
control parameters. Quantum states generated by processes like resonance fluorescence have low signal- 
to- noise ratio with poor control over the emission statistics [17]. Quantum optical processes like coherent 
superposition operations have been used to generate and characterize nonclassicality and entangled states 
[19]. The use of generic quantum optical models, such as semiclassical three-level atomic systems, that can 
be experimentally implemented and observed, can serve as an important tool to generate, investigate, verify 
and control nonclassical correlations and their features. 

In this paper, we investigate the nonclassical correlation properties of the photon states emitted from a 
three-level atomic system interacting with two classical driving fields. The interactions generate two-mode 
single photon states, arising from two controlled coherent transitions connecting the three levels, under the 
single photon approximation (SPA) [20]. The system can be set up in three different configurations, 5, A and 
V. We exhaustively study the correlation properties of the emitted two- mode photons, and characterize the 
dynamics of the different measures of nonclassical correlations. We establish a qualitative relation between 
the two different theoretical classes of correlation measures, namely, entanglement and the measurement- 
based correlations such as MID, QD and WD. The control parameters in the system enable us to define 
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specific regimes where certain correlations are enhanced based on the nature of the output photon states. We 
also analyze certain interesting features of the correlations generated by the interaction that throw light into 
the monotonicity and hierarchy of the set of measures used. The arrangement of the paper is as follows. We 
briefly discuss the different configurations of the three-level atom in Sec. 2. Then we have a short segment, 
in Sec. 3, defining the different correlation measures. In Sec. 4, we define the theoretical model used and 
the working approximations considered in the analysis. In Sec. 5, the numerical results obtained from the 
theoretical model are analyzed. We conclude in Sec. 6, with a summary of the results obtained and their 
possible ramifications. 

2. The three-level atom 

In this section, we briefly review our system. A three-level atom can be used in three different config- 
urations, namely, H, A and V [21, 22]. As a specific example, we focus on a gas of rubidium (Rb) atoms 
[23]. The energy levels 551/2, 5P 3 / 2 and 5£> 5 / 2 of Rb can be suitably used to generate each of the three 
configurations shown in Fig. l(a)-(c), as elaborated in the subsections below. Level 5Si/ 2 is the ground state 
and does not decay. Level 5D 5 / 2 is metastable, with a decay rate of about 1.0 MHz. The decay of levels 
5D 5 / 2 and 5P 3 / 2 are at rates of s» 1.0 MHz and Tp ps 6.0 MHz, respectively [23]. The conditions on the 
driving field Rabi frequencies f2j shown in Fig. l(a)-(c) are explained later in Sec. 4.3. 

2.1. The H system 

The cascade 5 system (Fig. 1(a)) uses the allowed dipole transitions |1) «-» |2) and |2) f-> |3), with two 
classical fields of Rabi frequencies Q\ and Q 2 driving these transitions, respectively. The field detunings 
are Ai and A 2 , set to zero throughout our analysis for near-resonant fields. The transition |1) f-> |3) is 
dipole forbidden. Levels |1), |2) and |3) correspond to the atomic levels 5<Si/ 2 , 5P 3 / 2 and 5P 5 / 2 of the Rb 
atom, respectively. Thus the decay rates of |3) and |2) are T 2 = Tp, = 1.0 MHz and Li = Tp = 6.0 MHz, 
respectively. Level P 3 / 2 serves as the shared level |2) during the interaction. The initial atomic state is 
ground state (|1)) populated and the levels |2) and |3) are unpopulated. S systems have been extensively 
used in coherent population trapping [24] and also in experiments to achieve laser cooling in trapped ions 
[25]. 

2.2. The A system 

The A system configuration can be obtained by folding the S, with levels 5S'i/ 2 , 5P 3 / 2 and 5D 5 / 2 of the 
Rb atom now marked as levels |1), |3) and |2), respectively, as shown in Fig. 1(b). With this identification for 
Rb, we observe that level |2) is energetically higher than level |3). This corresponds to a negative transition 
frequency u> 2 . The rotating wave approximation (RWA) thus holds for the negative frequency term of the 
field in the Hamiltonian [23]. Hence, the transition from level |3) to |2) annihilates a photon instead of 
creating a photon. The allowed dipole transitions are now |1) «-» |3) and |2) f-> |3), with two driving fields 
with Rabi frequencies fii and f2 2 now acting on these transitions. Level |3) is the shared level, and the 
transition |1) <H> |2) is now dipole forbidden. The decay rates of |2) and |3) are Ti = Tp, = 1.0 MHz and 
r 2 = Tp = 6.0 MHz, respectively. The initial atomic system is again ground state (|1)) populated, and 
the detunings are taken to be zero. The interactions of the three levels are distinctly different from the S 
system, and hence can be associated with different nonclassical behaviors. A systems have been extensively 
used in demonstrating diverse coherent phenomena such as stimulated raman adiabatic passage [26] and 
electromagnetically induced transparency (EIT) [14]. 

2.3. The V system 

The configuration of the V system (Fig. 1(c)) is considerably different from the 5 and the A systems. 
This is due to the fact that the shared level in the V system is the ground state. For the V system using Rb, 
we consider level 5Si/ 2 as the shared ground level (|1)) and two hyperfine levels of 5P 3 / 2 as the two-excited 
levels (|2) and |3)). Hence, T x = T 2 =T P = 6.0 MHz. The allowed transitions are |1) O |3) and |1) o |2), 
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driven by the classical fields of Rabi frequencies fii and 0,%, respectively. The transition |2) O |3) is dipole- 
forbidden, i.e., the ground state excitations take the system to two excited levels that cannot be coupled, 
and interactions are thus limited to ground state transitions. The initial atomic system is again ground 
state (|1)) populated, and the detunings are set to zero. V systems are widely used to study nonclassical 
phenomena such as quantum jumps [27], quantum Zeno effect [28] and quantum beats [16]. 

3. Quantum correlations 

We compare measures of quantum correlations that are defined from two different perspectives. Information- 
theoretic measures, such as MID [10], QD [11, 12] and WD [13], are based on the modification of quantum 
correlations upon measurement. These correlations calculate the difference in some specific property be- 
tween quantum states and their measured classical projections to give us a measure of the nonclassicality. 
On the other hand, there are nonclassical measures based on the entanglement-separability criterion. We 
use concurrence [29] as an example of entanglement monotone [8] . 

3.1. Measurement based correlations 

Measurement Induced Disturbance (MID): It is derived from the understanding that a truly classical state, 
with respect to some measurement, will remain unchanged after the measurement [10]. Let us consider a 
bipartite density matrix p a b- If Bf and Bj are complete von Neumann measurements (one dimensional 
projections) for subsystems a and b, respectively, for a classical state, 

Pab =J2 B ?® B ) P ab B i ® B r W 

ij 

The states p a b that do not satisfy (1) are essentially quantum in nature. MID measures the quantumncss 
in a bipartite state p a b by measuring the difference in the quantum mutual information between the state 
p a h and its least disturbed classical state obtained by the measurement, p c iass = J2i &i Pab where Bi are 
the spectral projections of the state p a b- Thus 

M(p ab ) = 2(pab) - Z(pclass), (2) 

where I(p a b) = S(p a ) + S(pb) — S(p a b) is the quantum mutual information [30]. S(p) — —tr(plog 2 p) is the 
von Neumann entropy of a quantum state p. p a and pb are the reduced density matrices of the subsystem 
a and b, respectively. A4(p a b) is the nonclassical measure of MID [10]. Unlike other measures of correlation 
based on projective measurements, MID does not introduce any optimization on the measured states. Hence, 
MID serves as an upper bound on other measurement based correlations [31]. 

Quantum Discord (QD): It is defined as the difference between two classically equivalent expressions 
for mutual information when extended to the quantum regime [11, 12]. For the density operator p a b, the 
expressions for quantum mutual information are 

I{Pab) = S(p a ) + S{p b ) - S{p ab ) 

+ J(pab) = S(pa) - S(p a \ b ), (3) 

where S{p a b), as defined earlier, is the von Neumann entropy. I(p a b) is the quantum mutual information [30] 
and S(p a \b) is the quantum conditional entropy [32]. The two expressions in (3) are equal in the classical 
regime. 

To calculate the quantum conditional entropy, S(ft a \b), we use the one-qubit orthonormal projection basis: 
= cos | |0) + sin § |1) and \i 2 ) = e -i * sin f |0) + cos § |1), with (i\j) = %, i,j = 0, 1. The projection 
density operators are given by the relations B\ = |ii)(£i| and B 2 = |ia)(*2|- On measuring the subsystem b 
using the above set of projection operators, the post-measurement states are p ab = — (I a <g> Bi p a b I a <8> Bi), 
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where pi = tr f, (I a C3> Bi p ab E Q (g)_Bj), and E Q is the identity operator acting on the subsystem a. The quantum 
conditional entropy is then given by 

S(p a \b) = ( 4 ) 

QD can then be defined using relation (3) and (4). 

J{Pab) = S{p a )-mm^piS{p l ab ), (5) 

i 

QD( Pab ) = I(p ab )-J(Pab)- (6) 

Quantum work deficit (WD): This concept is based on the fact that information is a thermodynamic 
resource [13]. It is defined as the amount of work (in terms of pure states) that can be extracted from 
a quantum bipartite system under a closed global operation and the amount that can be extracted using 
closed local operations and classical communications (CLOCC) [33]. WD is the deficit in the two operations 
due to loss of nonclassical correlations while performing CLOCC. Hence, WD is a measure of nonclassical 
correlations. 

Using the projection basis and the expressions for entropy defined previously, we can obtain the ex- 
pressions for the global and CLOCC operations. Under the class of global operations, the amount of work 
extractable (in terms of pure states) is given by la(Pab) = log 2 dim'H — S(p ab ), where AimH is the dimension 
of the Hilbert space. The amount of work that can be extracted using CLOCC is dependent upon local 
unitary operations, local dephasing and classical communication of the dephased state. p a b —> Pab Bi. 

Pab = Y.( Ia B * P ab Ia ( ? ) 

i 

The amount of work that can be extracted using CLOCC is given by iL(pab) = log 2 dirnH— infciocc^CPab)]- 
WD is then defined by the following expression: 

A( Pab )=I G ( Pah )-I L (p ab ). (8) 



3.2. Entanglement measure 

Concurrence: Entanglement in a bipartite system can be measured using this entanglement monotone 
defined for mixed states of two qubits [29] . Concurrence can be defined for a two-qubit density matrix p{t) 
as C(p) = max[0, Xi — A 2 — A 3 — A 4 ], where \(i — 1,2,3,4) are the square roots of the eigenvalues of the 
spin-flip operator, R = p(t)p(t), with p(t) = (a y £g> a v )p(t)(a v ® a y ), and a y is the Pauli spin matrix, 




Concurrence is closely related to the entanglement of formation. The entanglement of formation is a mono- 
tonically increasing function of concurrence. Concurrence, however, is not a resource-based measure such as 
entanglement of formation [34]. 



4. Model 

Three-level atomic systems have been extensively used to study quantum and nonlinear features of the 
semiclassical atom-field system [21] and also the nonclassical nature of emitted radiation in such systems 
[35] . Entanglement properties of three- level atomic systems have also been investigated [36] . However, a 
general classification of information-theoretic correlations for such systems does not exist, and we wish to 
investigate and compare the features of the above quantum correlation features using this versatile system. 
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4-1. The Hamiltonian 

The Hamiltonian for a general three-level atom interacting with two classical driving fields, in the RWA, 
can be written as [16]: 

H = Ho+Hi, (9) 

-Ho = fiw U |l)<l|+fi«22|2)<2|+fiW33|3)<3|, 

Hi = -h/2 (^ 1 e- l ^e- l ^ t \m)(n\+Vl 2 e- l ^e~ 1 ^ t \l){k\+ 1 R.c.), (10) 

where fax) a is the energy of level \i) (i = 1,2,3); fX,-e _l< ^ (j = 1,2) is the complex Rabi frequency corre- 
sponding to the classical driving field of frequency Vj, to, n,l,k = {1, 2, 3} denote the three atomic levels as 
appropriate for the S, A or V configuration. For the S configuration, (to, n, I, k) = (2, 1, 3, 2) correspond to 
the atomic transitions |1) o |2) and |2) O |3). For the A configuration, (m,n,l,k) — (3,1,3,2) correspond 
to the transitions |1) o |3) and |2) «-» |3), and for the V configuration, (m,n,l,k) — (3, 1,2, 1) correspond 
to the transitions |1) O |3) and |1) |2) (Fig. l(a)-(c)). 

4-2. The Atomic Density Matrix 

The state of the atomic system, at any time t, can be written in the following form: 

\ip(t)) A = C^e-^+^ll) + C 2 (t)e- l(w22+ « 2)t |2) + C 3 (t)e- l(w33+ « 3) *|3), (11) 

where = 1, 2, 3) are phases that depend on the detunings in a specific configuration. The detunings are 
defined as 

Ai = vi-ui, 

A 2 = V2 — W2, 

where oj% = u> mm — u> nn , and W2 = u>u — LUkk, are the transition frequencies, (m, n,l,k) have been defined 
earlier and are different for the three configurations. 

Using the wavefunction (11), one can create a pure state atomic density matrix pA(t) which depends on 
the classical driving field frequencies Vj{j — 1,2). We take phenomenological parameters to denote spon- 
taneous decays of the excited atomic levels [16] . The parametrized decay terms may lead to spontaneously 
generated coherences in the decay paths [37]. Such phenomenological decay models are commonly used in 
the study of quantum features in EIT [23], in quantum state tomography of emitted field states [17], and in 
earlier studies on atom-photon entanglement [38] . With the level decay terms, the dynamics of the system 
is in general mixed, and can be obtained using the von Neumann (quantum Liouville) equation of motion, 

PA® = -~[H,p A (t)} - ~{T,p A (t)}, (12) 

where the elements of the relaxation matrix T are the decay rates, (i|r|j) = T^iSij, i,j — 1,2,3. The 
time-evolved mixed atomic density matrix can be obtained provided the initial states of the atom (before 
interaction) are known. 

4-3. The Photon Density Matrix 

The nonclassical nature of the emitted radiation is dependent on the interaction between the three-level 
atomic system and the two-mode classical driving fields. The desired output, in our case, is to limit the 
generation to single photons for the two modes emitted after the interaction, so that at any given time 
within the lifetime of the atom, there will exist two photon states for each mode. Thus the two-photon 
density matrix can be written as 

Pvh = ^2Pij,i'j'\ i j)( i 'f\> 

where \i,j) , j')) stands for the two-photon states, with i and j (i' and j') = 0, 1 representing the number 
of photons in the first and second modes, respectively. Such a two-photon state can be achieved using the 
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single photon approximation (SPA) [20] within the rotating wave approximation (RWA) [16]. The RWA 
ensures that a photon is created only when an atomic de-excitation takes place. The SPA is applied by 
ensuring that the excitation time (due to the driving field strength) is larger than the decay time. If SI is the 
driving field Rabi frequency and T is the atomic decay rate, we require 1/S1 > 1/r for the SPA to be valid, 
so that the time taken for an atom to excite is much greater than the decay time and for small times only 
a single de-excitation will occur generating a single photon. 1 Hence, the ground state excitation strength fl 
should be smaller than the upper-level decay rate T. 

The output state is thus a two-qubit (bipartite) photon state. For a semiclassical interaction involving 
atoms and driving fields, it has been shown that the density matrix of the output radiation state can be 
completely derived from the atomic density matrix [17]. Under the far-field approximation [16, 39], for an 
atom located at fo, the field operators of the emitted radiation at the point of detection fare proportional 
to the atomic spin operators at the retarded time (t — \ f— fo\/c). This equivalence leads to an expression for 
the photon density matrix p p h(t) that is identical to the atomic density matrix at an earlier time, pA(t — r/c), 
calculated using the von Neumann equation of motion. The photon density matrix has a reduced rank three 
[40]. It has been shown [17] that such an equivalence leads to the complete determination of the output 
photon state using quantum state tomography, where measurements can be made on either the atomic or 
the photonic operators. Since the photon states can be completely determined by the atomic density matrix, 
the coherence of the photon correlations is closely related to the evolution of the atomic state. 

5. Results and Analysis 

As mentioned earlier, the system we consider is a gas of Rb atoms. The three levels 55 l 1 / 2 , 5P 3 / 2 and 
5D 5 / 2 of the Rb atom are appropriated to obtain the S, A and V configurations. All rates and frequencies 
are rendered dimensionless by scaling with the metastable level decay rate (~ 1 MHz). The scaled decay 
rates of 5P 3 / 2 and 5D 5 / 2 are Tp = 6.0 and = 1.0, respectively, and 5Si/ 2 is the ground state (Ts = 0) 
[23] . For desired results of two- mode single photon generation, we restrict ourselves to regimes that satisfy 
the SPA. In the S and A configuration (Fig. 1(a) and (b)), since the shared level is 5P 3 / 2 , we take the 
ground state excitation field (Jli) to be always less than the decay constant of 5P 3 / 2 (in 5, fii < Fi = T P 
= 6.0; in A, Qi < T 2 = Tp = 6.0). For the V configuration (Fig. 1(c)), the ground state excitation leads 
to transitions to the hyperfine levels of 5P 3 / 2 , and hence both the driving fields fli and Jl 2 are less than 
the decay constant of 5P 3 / 2 (£li,2 < Ti i2 = Tp = 6.0). In our analysis, we set the atom-field detunings 
to zero [41]. We calculate the measurement-based quantum correlations and entanglement of the output 
photon density matrix that can be derived using the atomic density matrix (Sec. 4.1) obtained from the von 
Neumann equation of motion (12), in the three configurations. 

We discuss the interesting quantum correlations generated in the two-mode photon state in specific 
parameter regimes using the 5 configuration as the reference. Some general observations can be made that 
are consistent with known results: MID always serves as an upper bound on the other measurement based 
correlations, such as, QD and WD [31]. The evolution of MID with respect to concurrence can be varied 
using the control parameters. There are two specific parameter regimes that correspond to two different 
hierarchies in the photon correlation. In Fig. 2, we consider the 2 configuration in two specific regimes of the 
driving classical fields. The decay constants for the S configuration are Ti = 6.0, T 2 = 1.0. Hence, the driving 
field strengths are in the range (f2i,f2 2 ) < 6.0. The detunings and the Rabi frequency phases have been 
set to zero. In Fig. 2(a), we consider the regime where fii > Jl 2 . Fig. 2(b), corresponds to the field regime 
fii < £! 2 (Sli = 2.0, Jl 2 = 5.0). For > £1 2 , MID forms a non-monotonic upper bound on concurrence 
at times t > 1.0. For the field regime fti < £! 2 (ili = 2.0, 51 2 = 5.0), concurrence forms a monotonic 
upper bound on the measurement-based correlations. We observe that the behavior of the correlations is 
closely related to the dynamics of the populations (inset of Fig. 2) . The non-monotonic behavior of MID is 
associated with the population difference in the two photon modes 1 00) and |11). It is clear from the plots 



1 We consider evolution times much smaller than the lifetime of the atom, since the condition of single photon emission may 
not remain valid due to spontaneous decay over a longer period. 
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Figure 2: (Color online) The time evolution for correlation measures MID (red continuous), discord (blue circles) and work 
deficit (green squares) along with the entanglement measure concurrence (black dashed) for the cascade (E) configuration. The 
field detunings are Ai = A2 = 0, and the phases of the Rabi frequencies are <f>i = <f>2 = 0. The level decay rates are Ti = 6.0, 
T2 = 1.0. SPA for this configuration requires that Q± < T±. The driving field strengths are (a) Qi = 2.0, CI2 = 1.0, and (b) 
f2i = 2.0, CI2 = 5.0. The inset shows the evolution of population elements of the two-photon density matrix and its purity. 



that the sudden increase in MID occurs when the populations of the modes 1 00) and |11) are nearly equal. 
This could be due to the fact that the non-optimization of the correlation measure in MID is skewed in these 
regions. For cases where MID is monotonic with the other measurement based measures, the population is 
distinctly unequal. Observing the purity in these regimes, one can state that the monotonicity is observed 
at higher levels of purity. 

A similar dichotomy in behavior can also be observed for fixed time dynamics of the system if the 
interaction is allowed to vary across driving field strengths. In Fig. 3, keeping the evolution time fixed (t = 
1.0) and varying the two classical driving field strengths, a similar behavior of the correlations is observed. 
MID is greater than concurrence and non-monotonic at times where the population levels are equal with 
significantly lower purity (Fig. 3(a)) as compared to regimes with unequal populations and higher purity 
where the measurement based correlations are monotonic to concurrence (Fig. 3(b)). Hence, we observe 
that the fixed time dynamics allows us to manipulate the correlation hierarchy by changing the ground-state 
driving field strength, Q%. The generation of monotonic correlations can be controlled by using parameter 
regions that allow higher purity in the output photon state. 

Similar parameter regimes can also be generated in the A and V configuration as shown in Fig. 4 and 
Fig. 5, respectively. Interestingly, the monotonic nature of the correlations in the A and V configuration is 
different from that of the H configuration. In the relatively high ground state excitation regime (high ill) 
in Fig. 4 we observe that the correlations attain steady-state values faster than in the S configuration. The 
measurement-based correlations such as QD and WD are not monotonic with either concurrence or MID 
at smaller times unlike in the S configuration where concurrence is always monotonic with QD and WD. 
There is a temporal discontinuity of concurrence around t &s 0.5. The concurrence collapses to a small finite 
value before reviving sharply. The revival of entanglement is associated with an increase in discord in the 
vicinity of the collapse. Such a feature of the correlations has been reported in other systems [42]. Other 
measures do not exhibit any discontinuity. At greater times (t > 1.0), the correlations are steady and weakly 
monotonic. The behavior of the correlations is again related to the population dynamics and purity of the 
density matrix as evident from Fig. 4 and Fig. 5. 

The different behavior of the correlation monotonicity in the A and V configuration as compared to 5 can 
be understood from the structural difference in the arrangement of the atomic levels. The highest excited 
level in the S configuration is the metastable state with a decay rate T2 ~ 1.0. In contrast, the A and V 
configurations have T2 ~ 6.0. Hence, the evolution of population dynamics and the temporal steady state 
occurs faster (t ss 0.5) than in the 2 configuration. However, the steady state bounds of MID or concurrence 
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Figure 3: (Color online) Fixed time (t = 1.0) MID (red continuous), discord (blue circles), work deficit (green squares), and 
concurrence (black dashed) in the S configuration as a function of the driving field strength The field dctunings are Ai 
= A2 = 0, and the phases of the Rabi frequencies are (f>i = 4>2 = 0. The level decay rates are Ti = 6.0, T2 = 1.0. SPA for 
this configuration requires that Q± < Ti. One driving field strength Qi is fixed at (a) 1.5, and (b) 3.5. The inset shows the 
variation of population elements of the two-photon density matrix and its purity. 



in different parameter regimes are common to ali the configurations. From Fig. 5, we observe that the steady 
state population dynamics in the V system results in high purity. We get an overlap of all measurement 
based correlations monotonically bounded by a low concurrence. This is due to the uniform decay rates of 
the two excited levels leading to a uniform distribution of population probabilities between 1 00) and |11) 
with negligible population in |01). 

In Fig. 6, we present a comparative study of the behavior in the three configurations under investigation, 
by choosing a common parameter regime that satisfies the SPA for all configurations. We consider a 
region of moderately low values of the driving fields (Hi = ^2 = 2.0) and another region of higher values 
(f2i = il 2 = 4.0). Some of the aspects of the correlation that can be qualitatively studied are monotonicity, 
temporal steady state, qualitative hierarchy and the nature of the two-photon density matrix. 

The correlation dynamics of the S system is dominated by the population dynamics of the photon state 
which evolves relatively slow due to the metastable nature of the highest excited level. We observe that the 
correlations do not achieve temporal steadiness in the observed times. This also leads to a lack of genuine 
monotonicity, with MID forming an upper bound in both the high and low field regimes. In comparison, 
the A and V systems have less stable excited states and hence achieve steady state correlations faster. For 
the A system, steady state correlations are fairly monotonic with concurrence bound at lower fields and 
MID bound at higher fields. V systems have relatively low values of steady-state quantum correlations, as 
compared to the long time values in S systems, with a concurrence bound at all fields. 



6. Summary 

The generation of mixed state quantum correlations is an important problem in the context of future 
implementation of quantum information tasks. The characterization of such correlations in a possible exper- 
imentally realizable model is thus an worthwhile exercise. With developments in experimental techniques 
and measurements of information-theoretic correlations like QD [43] , the importance of generating and ana- 
lyzing correlations has manifestedly increased in recent times. The possibility of using well studied quantum 
optical systems like the three-level atom to generate, characterize and parametrically control mixed state 
quantum correlations is indeed an encouraging step in this direction. 

In this paper, we have used a semiclassical three-level atom interacting with classical driving fields to 
generate a correlated two-mode photon pair controlled by the driving field parameters and the phenomeno- 
logical decay terms of the atomic level. Under certain physical approximations, we have exhaustively studied 
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Figure 4: (Color online) The time evolution for MID (red continuous), discord (blue circles), work deficit (green squares) 
and concurrence (black dashed) for the A configuration. The field dctunings are Ai = A2 = 0, and the phases of the Rabi 
frequencies are <f>i = (p2 = 0. The level decay rates are Ti = 1.0, T2 = 6.0. SPA for this configuration requires that f2i < T2. 
The driving field strengths are thus taken as (a) Qi = 4.0, Q2 = 5.0, and (b) f2i = 4.0, Q2 = 1-0. The inset shows the evolution 
of population elements of the two-photon density matrix and its purity. 




'0 05 i L5 2 15 3 
t 



Figure 5: (Color online) The time evolution for MID (red continuous), discord (blue circles), work deficit (green squares) 
and concurrence (black dashed) for the V configuration. The field dctunings are Ai = A2 = 0, and the phases of the Rabi 
frequencies are ipi = <l>2 = 0. The level decay rates are Ti = T2 = 6.0. SPA for this configuration requires that fix 2 < Pi, 2- 
The driving field strengths are Qi = 2.0, Q2 = 4.0. The inset shows the evolution of population elements of the two-photon 
density matrix and its purity. 
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Figure 6: (Color online) Time evolutions of MID (red continuous), discord (blue circles), work deficit (green squares) and 
concurrence (black dashed) for all the three configurations, S (top), A (middle) and V (bottom). The field detunings are Ai 
= A2 = 0, and the phases of the Rabi frequencies are <j>i = <p2 = 0. The chosen driving field strengths of f2i = Q2 = 2.0 (left 
panel) and Qi = Q2 = 4.0 (right panel) satisfy the SPA for all three configurations. 



and compared the dynamical behavior of various nonclassical correlations of the two-photon state generated 
in the semiclassical atom+field system. We have probed the dynamics of entanglement (concurrence) as well 
as of measurement-based correlations, such as MID, QD and WD, for three different configurations, namely, 
S, A and V, of the three-level atom driven by two controlled external classical driving fields. The qualitative 
characterization of the correlations based on the monotonicity, general hierarchy and steady state behavior 
are achieved using the field parameters and the decay terms. 

The correlation behavior is observed to be configuration dependent. The manipulation of control pa- 
rameters in different configurations leads to variations in the dynamic evolution of the correlations. The 
S configuration produces photon states with relatively high correlation even at low driving fields. The dy- 
namics in the S system is more mixed as compared to the A and V systems. The 5 system can be suitably 
controlled using the driving fields to generate correlations dominated both by MID or concurrence and is 
ideally suited to experimentally study the temporal evolution of the two measures with respect to the evolu- 
tion of the system in the Hilbert space. The optimization of MID and monotonicity of measurement based 
correlations can be experimentally analyzed using quantum state tomography. A and V systems are better 
suited for generating steady monotonic correlations in both low and high strength driving field regimes. 
The A system can be suitably tuned to generate steady correlations with either entanglement or MID as an 
upper bound. The mixedness in the generated states can be controlled using the driving field strengths for 
implementation in experiments. V systems, however, can generate ideally pure correlated photons bounded 
by concurrence at all field strength regimes. The absence of a metastable state in the V system allows 
production of pure correlated output photon states. The measurement-based correlations are all equal at 
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steady values. However, in the A and V systems, significant correlation is generated only for high driving 
field strengths. 

Hence, specific regimes and configurations can be used to generate and manipulate the correlations in the 
output two-photon state as desired. Our findings may thus be immensely useful in practical implementations 
with such interacting photon states. 
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